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0\ ■ Abstract 

On 

Nuclear processes involving momenta much below the mass of the pion may 
be described by an effective field theory in which the pions do not appear 
I/-) ■ as explicit degrees of freedom. The effects of the pion and all other virtual 

hadrons are reproduced by the coefficients of gauge-invariant local operators 
■^i- . involving the nucleon field. Nucleon-nucleon scattering phase shift data con- 

strains many of the coefficients that appear in the effective Lagrangean but 
at some order in the expansion coefficients enter that must be constrained by 
^^ ■ other observables. We compute several observables in the two- nucleon sector 

£— : ■ up to next-to-next-to leading order in the effective field theory without pions, 

Q^ . or to the order at which a counterterm involving four-nucleon field operators 

is encountered. Effective range theory is recovered from the effective field the- 
ory up to the order where relativistic corrections enter or where four-nucleon- 
external current local operators arise. For the deuteron magnetic moment, 
quadrupole moment and the np — ► d^y radiative capture cross section a four- 



in 



o 

nucleon-one-photon counterterm exists at next-to-leading order. The electric 
polarizability and electric charge form factor of the deuteron are determined 
up to next-to-next-to-leading order, which includes the first appearance of 



X 



$h ' relativistic corrections. 
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I. INTRODUCTION 

Weinberg's pioneering work [I|] on the application of effective field theory (EFT) to 
nuclear physics has given rise to a decade long effort to construct a perturbative theory of 
nuclear interactions @- f3lf . The small kinetic and potential energies of nuclear systems 
compared to the typical scale of strong interactions, either Aqcd or the chiral symmetry 
breaking scale A x , suggests that the ratio of these two scales might provide a small expansion 
parameter. Indeed, in the single nucleon sector the ratio of external momenta and quark 
masses, m q , to A x allows for a systematic field theory description of pionic and external 
current interactions with a single nucleon. However, extending this simple idea to multi- 
nucleon systems has proved to be a challenge. 

It is important to understand nuclear physics directly from quantum chromodynamics 
(QCD) and constructing the correct low-energy effective field theory is a central component 
of such a program. The optimal way to determine the properties of multi-nucleon systems is 
to match onto the low-energy effective field theory and perform computations directly with 
the effective theory, and not with the theory written in terms of quark and gluon fields. 
Further, one wants to be able to make such a matching between QCD and the low-energy 
effective theory free of any unjustified assumptions about the dynamics or structure of the 
hadrons. Unfortunately, at this point in time such a matching is not possible. However, 
the symmetries of QCD along with the fact that we wish to examine observables involving 
small external momenta allow the low-energy effective field theory to relate the amplitudes 
of different strong interaction processes, at some level of precision. 

Significant progress in understanding both the two-body and three-body systems in terms 
of effective field theory has been made during the last year. Several observables in the two- 
nucleon sector have been explored in detail in the theory with dynamical pions, such as 
the electromagnetic form factors of the deuteron |I2[ , jd — > jd [[L3|,|n]], np — > dj |T6| , p!7|j2^] 
and vd — > ud, vnp |K| . In addition, the inclusion of coulomb interactions into the pionless 



theory has been shown to be straightforward p6|. In the three- nucleon system it has been 



demonstrated that low-momentum processes can be accurately and precisely described in 
the effective field theory without dynamical pions [|3l|. For example the quartet N — d 



scattering length has a rapidly converging expansion in terms of the ratio r$ jo^ Sl \ where 
7q is the effective range and a 1 - Sl > is the scattering length for A^ A" scattering in the 
spin-triplet channel. Calculations in the two-nucleon sector [[TZ|-|iT| with dynamical pions 
using the KSW power counting scheme [|I0j, (where pions and higher dimension four- nucleon 
operators are treated in perturbation theory) have been performed at leading order (LO) 
and next-to-leading order (NLO), and just recently been extended in some cases to next-to- 
next-to-leading (NNLO) 1J|L9|]. 

In this work we make a thorough investigation of the two-nucleon sector in the theory 
without dynamical pions, appropriate for very low-momentum processes. Effective range 
theory [^2J is found to reproduce the effective field theory up to the order where relativistic 
corrections or local four-nucleon-external-current local operators enter. In many cases effec- 
tive field theories for hadronic systems are used to systematically include the approximate 
SU(2)i ® SU(2)r chiral symmetry of QCD. As we are considering the theory of nucleon 
interactions where pions do not explicitly appear, chiral symmetry does not play a central 
role and the usefulness of this investigation could be questioned. One motivation for con- 



structing and understanding this theory is that calculations of the observables considered 
in this work with any scheme, such as potential models or the effective field theory with 
dynamical pions, must reproduce the analytic results we obtain. A secondary motivation is 
that this theory provides a very clear demonstration of how the perturbative effective field 
theory approach to NN interactions behaves at higher orders, including relativistic effects. 

II. NN SCATTERING IN THE 3 S"i - 3 £>i CHANNEL 

The large scattering lengths in the two-nucleon system require a modification of the 
power counting rules that follow naturally from the single nucleon sector in order to de- 
scribe multi-nucleon processes. For processes in which all external momenta are much below 
the mass of the pion (|k| <C m n ) it is desirable to use an EFT without pions, which we 
denote by EFT(tF). The external momenta divided by the pion mass is the small expansion 
parameter of this theory, denoted by Q. As there are no dynamical pions, all quark mass 
effects are reproduced by local operators involving two or more nucleons and an arbitrary 
number of spatial derivatives. EFT (ft has the same power counting as effective range theory 
where the relative momentum of the nucleons in the two-body NN system is treated as the 
small expansion parameter. However, EFT(/f) includes operators with insertions of arbitrary 
numbers of external currents, that are not constrained by NN scattering phase shifts alone. 

Let us recall the low-momentum behavior of the phase shift in a system with a large 
scattering length. We know from elementary scattering theory that |k| cot#o is an analytic 
function of external kinetic energy with a radius of convergence bounded by the threshold 
for pion exchange in the t-channel, and thus 

|k|cot(J = — + 2 r °l k ! 2 + ri l k ! 4 + - ' ( 21 ) 

where we have chosen to write the expansion in terms of the center of mass momentum 
|k|. One is free to expand |k| cot^o about any point within the region of convergence, and 
eq. (|2.f|) corresponds to expanding about |k| = 0. The size of the coefficients ro,ri, ... is 
set by the range of the interaction between the nucleons, of order the pion mass, and the 
Ti are taken to scale as Q° in the power counting. In contrast, the scattering length a is 
not constrained by the range of the underlying interaction. The scattering lengths for NN 
scattering in both the ^o and 3 S± channel are very much larger than the inverse pion mass, 
and are taken to scale as a ~ Q^ 1 in the power counting. In dealing with the deuteron 
bound state, it is convenient to make the expansion of |k| cot<5o about the location of the 
deuteron pole f32"f , and not about |k| = 0. Expanding around the deuteron pole, |k| 2 = —7^ 
gives 

|k|cot«J = -7t + 2^(l k | 2 + ^) + ^2 (|k| 2 + 7 4 2 ) 2 + - 

= 0(Q) + 0(Q 2 ) + 0(Q i ) + ... , (2.2) 

with 7 ( _1 = 4.318946 fm, pd = 1.764 fm, and W2 = 0.389 fm 3 [33]. The total energy of each 



nucleon in the center of mass is En = y |k| 2 + M N , and the deuteron pole is located at 
|k| = ij t , the solution to 



-B = 2y/M%-tf-2M N , (2.3) 

where B = 2.224575 MeV is the deuteron binding energy. Subsequently, we will find it 
convenient to work with the parameter 7 = y/M^B, which is related to 74, by a Q expansion 



of eq (|J), 



7 = * + jjjjr + - • P-4) 

The S-matrix describing scattering in the coupled channel J = 1 system is written as 

/ e i25o cos2ei ie i ^ So+52 hm2e 1 \ , , 

Ue i(<5o+<52) sin2e 1 e i252 cos2e 1 J ' ^ ' 



where we use the "barred" parameterization of |p4| , also used in ||33[ . It will follow naturally 
from the EFT(/F) that S\ is supressed by Q 2 compared with 6q , and therefore, up to N 4 LO, 
we can isolate the S-wave from the D-wave in the deuteron channel, leaving 

Soo = e m ° = 1 + — ^— . (2.6) 

cot o — I 

The phase shift 5 has an expansion in powers of Q, 5o = 5$ + 5q + 5q + ..., where the 
superscript denotes the order in the Q expansion. By forming the logarithm of both sides 
of eq. ( |2.6| ) and expanding in powers of Q, it is straightforward to obtain 

^ 0) (|k|)=7T- 




# ) (|k|) = -f|k| 



^ 2) (|k|) 



Pdl , 7 3 



+ 



4 8ikT^(7 2 + |k| 2 ) 



(2.7) 



which are shown in fig. (|XJ) . Up to NNLO, the shape parameter term, W2 does not contribute 
to the S-wave phase shift and therefore, in addition to being numerically small (about a 
factor of 5 smaller than pd), W2 enters only at high orders. Formally, the first deviations 
from linear |k| dependence at large momenta arise from relativistic corrections (the second 
term in Sq in eq. Q2.7]). The Q expansion is clearly demonstrated in eq. (|2.7|), by simply 
counting powers of |k| and 7, which both scale like Q. We expect that this perturbative 
expansion of the phase shifts converges up to momenta of order ~ m n /2 ~ 70 MeV, at which 
point one encounters the t-channel cut from potential pion exchange. In the theory where 
pions are dynamical, the range of convergence of the theory is much beyond the scale set by 
the pion mass. An analogous expression for the perturbatively expanded S-wave phase shift 
that is valid up to momenta of order ~ 300 MeV at NLO can be found in ||10|| . 

This simple perturbative expansion of the phase shift and scattering amplitude is repro- 
duced order by order in the EFT(^) expansion. Relativistic corrections are encountered at 
NNLO and therefore a discussion of the nucleon two-point function is appropriate. The total 
energy of a nucleon E N and its kinetic energy T N are simply related by E N = M N + T N . A 
nucleon propagator can be written in terms of T/v and the nucleon three momentum p, 
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FIG. 1. The phase shift 5q as a function of the center of mass momentum |k|. The dashed 

the solid curve corresponds to 



curve corresponds to 5q , the dotted curve corresponds to 5$ +5$ 



;(°) j_ *(1) 



?(2) 



Sq + 5q + 5q , and the dot-dashed curve is the Nijmegen partial wave analysis [St] . 
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(2.8) 



and is reproduced by a Lagrange density involving the nucleon field N with the time- 
dependent phase corresponding to the nucleon mass removed, 



C = iV f 



id + 



n 2 



2M, 



N 



2M 



N 



N + 



(2.9) 



In the EFT(/f) we treat the third term in eq. fl2.9|) as a perturbation and it is not resummed 
into the nucleon propagator. In order to recover the usual nonrelativistic expansion of 
the two-point function, a field redefinition is employed followed by use of the equations of 
motion p5|, e.g. 



N' 



1- 



4M^ 



+ ...IN 



(2.10) 



gives 



C = N 



't 



id + 



2M 



N 



W% 



+ 



N' + 



(2.11) 



It is somewhat inconvenient to use this second form of the Lagrangean due to the fact that 
the field redefinition in eq. (|2.10| ) must be performed on all terms in the Lagrange density 
describing the NN systems including the multi-nucleon interactions. 



The inclusion of electromagnetic interactions requires that V operators be replaced by 
covariant derivatives D = V — ieA, and time-derivatives do be replaced by D = do + ieA . 
The Lagrange density describing the interaction between two nucleons scattering in the 3 5i 
channel written in terms of one-body C\ and two-body £ 2 interactions, and up to NNLO, is 



D 2 



D 2 



d = iV f iD , 

[ 2M N 2M N 

C 2 = -fcf Sl) (iV^JV) 1 (N T PiN) 



N 



+ 



i C ( i Si)_ 



(N T PiN) ] [N 



PiT? + rJ 2 p - z\5Pit) 



16 



*ci 3Sl) (n t \pT) 2 + BV, - 2rJ pji 



N) IN 



tT 



N) +h.c. 

2 4—2 



prj +15 Pi- 2 i 5p i i) 



N) (2.12) 



where Pi is the spin-isospin projector for the 3 Si channel 



Pi = —7=r<J2<JiT~2 



Tr P}P, 



-Oi 



(2.13) 



The subscript on the coefficient denotes the number of derivatives in the operator. There is 
another operator involving four derivatives that we have not shown in eq. (|2.12| ), 



L ~ "32 ° 4 



(iV T (lD - tifPiN^ (N T PiN 



+ h.c. 



(2.14) 



(we have not shown the correct ordering of the Pi and D operators). Renormalization group 
scaling |10[ of the operators in eq. ( |2.12| ) and eq. (|2.14j ) indicates that while the contribution 



from 



•Sx) 



is NNLO, the contribution from nj\ is N 3 LO. The time-ordered product of 
two ^C*2 operators does not induce the momentum structure of the K2\ operator. 

The effective range expansion provides a complete description of scattering in the low 
energy region. It is straightforward to show that the relation between the expansion of cot 5q 
and the coefficients in eq. (|2.12|) is, 



Ikl cot<5o 



2nE 



M 2 N EtJ 2n \k\^ 



+ V 



(2.15) 



where fi is the dimensional regularization renormalization scale, explicitly introduced by the 
PDS subtraction procedure [|K] . 

E = 2M N + T 



2M, 



N 



M, 



N 



4M 3 , 



+ 



(2.16) 



is the center of mass total energy of the two nucleon system, where each nucleon has kinetic 
energy -|, and momentum of magnitude |k|. Lorentz invariance ensures that for two nucleons 
moving with total momentum P and kinetic energy T, these quantities are related to the 
center of mass quantities by 



T + 



21 

AM, 



N 



AM 



T 



N 



21 

AM, 



N 



(2.17) 



Notice that in eq. (|2.2|) momentum independent terms appear at each order in the Q 
expansion, as the momentum expansion is about the deuteron pole and not about k = 
(further, one sees contributions from all powers of |k| 2 , up to a maximum order dictated by 
Q). Consequently, for the effective field theory to reproduce the effective range expansion, 
the coefficients appearing in eq. ( |2.12| ) will have an analogous expansion in powers of Q, e.g. 

t r { 3 Si) _ ^( 3 5i) ^( 3 Si) ^( 3 5i) 
°0 _ °0,-l "T °0,0 "T °0,1 "T ••• 

^( 3 5i) _ ^( 3 Si) ^( 3 5i) 
( - y 2 ~~ ^2,-2 "T ^2,-1 i" ••• 

t7f 5l) = t^ + ... . (2.18) 

The second subscript on each coefficient denotes the powers of Q in the coefficient itself. 
Relating terms order by order in the k expansion of eq. ( |2.2|) and eq. ( |2.15|) we find that 
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M^ (// - 7) 
00,1 ~ M^-7) 3 + 2M JV M^- 7 )3 ' °-° 55tm 



0.559 fm 2 

2„,4 „ „,3 



^ 



2vr 

27T p 2 7 2 27T 1 



*® = #r^S = 10 - 420 fm4 

Mjv (/x - 7) 



ilM/i-7) 3 M N M^~l) 



-2.210 fm 4 



,2 



*® = -i^o^? = - 19 ' 440 fm ° • (2 ' 19) 

where the l/Mjf terms that do not depend upon pa are relativistic corrections. We have 
chosen to renormalize the theory at // = m^. Strictly speaking we should choose a much 
lower scale /j, ~ Q, but the RG invariance allows us to renormalize at any scale once we have 
established the power counting. The hierarchy between coefficients of the same operator 
but of different orders in the Q expansion is clear. It is straightforward to show that these 
coefficients reproduce the perturbatively expanded phase shift in eq. Q2.7|). 

The amplitude that one computes from Feynman diagrams is simply related to the S- 
matrix. Explicit computation of Feynman diagrams, in an arbitrary frame, gives 

_ 2n ^ l 

Anr ~ ~m% |k|cot<y-*|k| ' (2 ' 20) 

which is related to the S-matrix by 



2% |k|M 2 

i-NR 



S = l + "" . = 1 + i^^A NR . (2.21) 

cot 6 — i txE 

In the process of computing the S-matrix from the Feynman diagrams, we have used the 
following S-wave states (consistent with the sign convention of p4| ) 



|l,a;P,|k|) L= 



4vr (2tt) 3 



dQr, 



N% ,. P a N P 



—k 



f+k 



10) 



(2.22) 



for the J = 1, Jz = a and L = NN state. P is the momentum of the center of mass of the 
NN system. Similarly, for the D-wave states, 



|l,a;P,|k|) L= 



8tt |k| (2tt) 



dQi 



k x k a_i 5 ax, k |2 

3 ' ' . 



N% , P x Nr_ 



t 



10} 



(2.23) 



where the —ve sign gives the correct phase for the mixing parameter e%. These states are 
normalized such that 



L ,(l, b; P', |k'||l, a; P, |k|) L = <5 3 (P - P')<J(|k| - |k'|)^ 



LL' 



(2.24) 



Scattering between the S-wave and D-wave is induced by local operators first arising 
at Q 1 in the power counting. At order Q l and Q 2 , the lagrange density describing such 
interaction is 

4*0 = jt&**> (iV r P i JV) t (N T <D^' xyj N)T ijxy + h.c. 



16 



j c (sd) N 



tT 



PiT$ 2 + £> 2 fi - 2*5 pj5 



n) (N T O isd) ' xyj N)T i3xy + h.c. , (2.25) 



where 



Q(sd),xyj = ffffpj + pJf) X Y) y _ ffpij? - ffP'T? 

1 



rj-ijxy 



xixfijy 



.gijxzy 



n 



(2.26) 



and where n is the number of spacetime dimensions. When electromagnetic interactions are 



ignored the operator in eq. ( |2.26| ) collapses to 

(sd),xyj _^ pj fa _^\ x fa _^ 



^2.27) 



which is explicitly Galilean invariant. As the S-D mixing parameter is being computed up 
to NLO relativistic corrections do not enter. There is another operator with four derivatives 
that contributes to S-D mixing, 



Bf> = __€« (n t p'n 



*'jV 



'^\2Q(sd),xyj _|_ Q{sd),xyj fj^\2 _ ^f) Q{sd),xyjj^ 



N T lJxy 



(along with its hermitian conjugate) but it is of order Q 3 and does not contribute at NLO. 

The coefficients that appear in eq. fl2.26| ) themselves have an expansion in powers of 

Q, e.g. $Cq = *CqJi + ^Cq + .... Performing a Q expansion on the mixing parameter 

— —(2) —(3) 

£i = e\ + e\ + ... it is straightforward to demonstrate that 
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FIG. 2. The S-D mixing parameter ~e\, in degrees, as a function of the center of mass momentum 

—(2) 

|k|. The dashed curve corresponds to e\ , and the dot-dashed curve is the Nijmegen partial wave 
analysis. 



Renormalization group invariance of this leading order contribution to e\ indicates that 
^Cq _i oc (// — 7) -1 . Renormalizing at the scale /1 = m ff , and comparing with the Nijmegen 

phase shift analysis ||35|| , we find KJqJ\ = —4.57 fm . A comparison between the Nijmegen 
phase shift analysis and our fit can be seen in fig. (H). 
At the next order, Q 3 , the contribution to e\ is 



^r(iki) 



3 4?r 
+ 



(^^ + (^-7)^)7: 



/'y 2 



(2.30) 



where we have expanded *U£ = ^2-2 + •••• The first term in eq. ( |2.30| ) has the same 
momentum dependence as e\ , the leading contribution in eq. ( 2.29|) . Requiring that the 



_(2) 

momentum structure of e\ does not appear at higher orders in the expansion constrains 
*±°tobe 



"0,0 



Further, as the second term in eq. ( |2.30|) is an observable, we obtain the RG equation 

% (f *^ + (/i " 7) ^ ) = ° • (2 ' 32) 

By considering /1 to be of order the short-distance matching scale it is clear that the combi- 
nation |pd ^0-1 + (A* — 7) ^2-2 appearing in eq. (|2.30| ) and eq. (|2.32| ) is of order Q° and 



not Q 1 , as naively guessed from the subscripts of ^Cq _{ and ^2-2- Therefore 



^ 3) (|k|) = , (2.33) 

and the first corrections to s\ arise at NNLOQ, i.e. e[ . 

Local operators that contribute to D-wave to D-wave scattering contain four or more 
powers of the external momenta. A RG analysis shows that the coefficient of the leading 
interaction scales like /x _1 , through the time-ordered product of two *Cq operators. There- 
fore, contributions to the scattering amplitude start at order Q 3 (N 4 LO, i.e. 5^ is the first 
non-zero contribution to the 82 phase shift), which is higher order than we are working. 



III. PROPERTIES OF THE DEUTERON 

In the previous section nucleon-nucleon scattering was used to determine the coefficients 
that appear in EFT(yf) up to NNLO for S-wave interactions and up to NLO for mixing 
between the S-wave and D-wave. This provides us with enough information to determine 
several properties of the deuteron up to NNLO in EFT(yf). 

To begin with it is informative to recover basic kinematic properties of the deuteron from 
the effective field theory formalism. The inclusion of the deuteron into EFT was discussed in 



detail in |L2] , but relativistic corrections were not included. In order to compute the matrix 
element of the electromagnetic current between deuteron states, one starts by defining the 
3-point function for two nucleons with kinetic energy T interacting with the electromagnetic 
current, transferring four-momentum q and producing two nucleons with kinetic energy T", 

Gf j (T,T',g)= / d A x d 4 y e + ^ T ' x °- p '^e-< Tyo - p - y \0\T[V i (x)J^ m (0)V j (yy]0) . (3.1) 

where T>i(x) = N T PiN(x). The S-matrix element for this interaction is found by LSZ 
reduction of this Green-function 

(p',j\JL\P^} = {^z) 2 [G-\T)G- 1 (T')Gf J (T,T',q)L T _,__, , (3.2) 

polc> pole 

where T po i e is the location of the deuteron pole for two nucleons moving with three- 
momentum p, 

rp 1 pole I P I I _|_ ' 

pole AM N ~ AM N ~ ~M^ AM% 

9 112 2112 114 

t -J- 1 l p i , t\p\ _ IpI 4 , ,00^ 



The expression for T po i e has been expanded in powers of M N because it is the nucleon mass 
that will naturally arise from Feynman diagrams computed with the effective field theory. 
An analogous expression holds for T' le . The ellipses in eq. (|3.3|) denote terms of order Q 6 



1 This was pointed out to us by S. Fleming, T. Mehen and I. Stewart, and can be found in their 
recent preprint 



10 



or higher. The two nucleon two-point function, defined in [12| can be generalized to include 
relativistic effects, 



G(T) = I d 4 x e 4(Tx °- px) (0|T[P i (a:)P](0)]0) 



6 HY 



iZ(T) 



Tpoic + «e 



(3.4) 



The wavefunction renormalization factors are straightforward to compute from a single 
nucleon-nucleon bubble, which is explicitly lorentz invariant (up to the order one com- 
putes). An important modification arising from the inclusion of relativistic effects is that 

ti -7T1-- This leads to wavefunction 



the two-point function is an explicit function of T + 
renormalization factors of 

-*- pole 

Z = Z{T po \ e 



4Mjv 4Mjv 



2Mjv 



dE(T)/dT 



(3.5) 



iwi; 



pole 



where S(T) is the irreducible bubble presented in |12|]. We find that 



4 Z 






, 2 2 7 2 (7/i-57) 

1 + IPd + 7 Pd + 



IMUv-l) 



+ 



(3.6) 



which is lorentz invariant, but depends on the renormalization scale p. 

Evaluating the matrix element of the Hamiltonian between deuteron states 



(p',j\H\p,i) = (p',j\ / d 3 x 
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M, 
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v 



2M, 



{2irYS d {p-p')S i 



(3.7) 



The interactions and relativistic contributions combine together in such a way to reproduce 
the deuteron energy-momentum relation at the order to which the calculation is being per- 
formed. This is a generic feature, which must be present for the theory to be consistent. In 
computing any observable, one finds that the deuteron mass is recovered order by order in 
perturbation theory. The issue of gauge invariance for the deuteron bound state is also sim- 
ply addressed. The lorentz-invariant lagrange density that describes the leading low-energy 
interactions of the deuteron has the form 



C = d\ 



iD + 



D 5 



D 2 



di + 



(3.8) 



2M d 2M d 

and so reproducing the correct energy-momentum relation for the deuteron order by order 
in the expansion forces the matrix elements to be gauge invariant order by order in the 
expansion. In fact, direct calculation including the leading relativistic corrections gives 



C = d\ 



i(d + ieA ) + (V - ieA) 



AM 



+ 



N 



SM 3 N J 



(d + ieA c 



1 



V4M, 



A' 



dj + 



. (3.9) 



Order by order one recovers the correct matrix elements of J£ m to reproduce the couplings 
induced by eq. ( |3.8|) and also the Thompson limit for photon-deuteron elastic scattering. 
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A. Electric Polarizability of the Deuteron 

The electric polarizability of the deuteron, C(eo, has been investigated thoroughly with 
potential models P8f - [(47]. To a very high precision one obtains c*eo = 0.6328 ± 0.0017 fm 3 



with those techniques. This is no surprise because the electric polarizability is dominated 
by the long range behavior of the deuteron wavefunction. If the model is tuned to reproduce 
this component of the wavefunction the predicted electric polarizability should be very close 
to nature. In effective range theory the polarizabilities are assumed to be dominated by the 
asymptotic S-wave component of the deuteron wave function, 



* (er) <'Wmt^>V ■ < 310 > 



which yields an electric polarizability of 

ER _ aM N 1 



a 



E0 — or, 4 

2 2 



32 7 4 1 - 1Pd 
aM N 



32 7 4 



i + 7Pd + rri + ... , (3.ii) 



where a is the electromagnetic fine structure constant. Numerically, a E ^ = 0.6338 fm , 
which is very close to the value obtained by the potential model calculations. 

With the EFT (ft) we have computed ol E q up to NNLO, including relativistic corrections. 
As with the phase shifts, it has a perturbative expansion in Q, a E0 = «£o + a Eo + a Eo +•••> 
and we find 



(-4) , (-3) . (-2) aM N 

a E0 + ®E0 + a E0 = -T^T 



,2^2 



1 + IP, + IPa + 3^ 



0.377 + 0.153 + 0.062 + 0.0006 
= 0.592 fm 3 . (3.12) 

Numerically, the relativistic corrections are very small, two orders of magnitude smaller 
than the NNLO corrections from the four-nucleon interactions. The value of ctEo shown in 
eq. (|3.12| ) is within ~ 5% of that computed with potential models and with effective range 
theory. Despite making a small contribution to the numerical value of a^o this calculation 
demonstrates that relativistic contributions can be computed easily with the EFT. The 
operators that mix S-wave and D-wave (corresponding to the D-wave component of the 
deuteron in potential model language) make contributions to a E . They do not contribute 
to the scalar polarizability, a E0 , but do contribute to the tensor polarizability, a E2 ■ Such 
operators will contribute to a E Q at higher orders in the expansion. 

The four-nucleon operators in the EFT(^f) reproduce the contributions from effective 
range theory in addition to the relativistic corrections that arise at NNLO. However, the 
situation is different at higher orders. a E0 will receive contributions from P-wave inter- 
action between nucleons (a two- derivative operator that is not renormalized by the large 
S-wave scattering length). a E0 receives contributions from the polarizability of the nucleons 
themselves. It is only a E0 and higher that receive contributions form the four-nucleon-two- 
photon interaction, a counterterm for the polarizability of the deuteron, that can only be 
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determined from two-nucleon and two-photon processes. Extrapolating the relative size of 
the contributions seen in eq. ( p. 12 ) we conclude that the counterterm contribution will be 
of order ~ 0.005, much larger than the relativistic corrections, but still only ~ 1%. 

In the EFT with pions included as a dynamical field, the electric polarizability at NLO 

isra 



aM 



&E0 



N 



32 7 4 
aM N 
32 7 4 
aM N 

32 7 4 



1 + C 2 {fi) 



M Nl ({i - 7 ) 2 g\M Nl ml{?,ml + 16m nl + 24 7 2 



2vr 



i + cm M -^-^ + AM,. 



1 + IPd + 



2n 

g\M Nl 10 7 2 
47r/ 2 3 m 2 



127r/ 2 (m* + 2 7 ) 4 
_ 8 7 16 7 2 
3 m„- 3 m 2 



AixP 



+ 



(3.13) 



where gu is the pion-nucleon axial coupling constant, and / is the pion decay constant. 
From the expansion of a^o i n powers of j/m^ that appears in eq. ( |3.13|) it is clear that one 
recovers the effective range expansion, but with an additional contribution from the pions 
that is higher order in the momentum expansion. In the EFT(^), 7 /m 7r ~ Q while in the 
theory with dynamical pions it is 7 /m 7r ~ Q°. The residual contribution on the third line 
of eq. ( |3.13| ) proportional to g\ is therefore N 3 LO in the momentum expansion appropriate 
to the theory without pions. 



B. Electric Charge Form-Factor of the Deuteron 

The electric charge form factor of the deuteron is a very well measured object over a 
wide range of momentum transfers. Potential models and effective range theory reproduce 
the data very well in the kinematic regions where they are applicable. (For a recent and 
very comprehensive review, see p9| .) 

A deuteron with four-momentum p M and polarization vector e M is described by the state 
|p, e), where the polarization vector satisfies p M e M = 0. An orthonormal basis of polarization 
vectors et satisfies 



P^i 






-5, 



ij j 



E 



1 { Ml 9 



(3.14) 



where Md is the deuteron mass. It is convenient to choose the basis polarization vectors 
so that in the deuteron rest frame ef = <5f . Deuteron states with these polarizations are 
denoted by |p, i) (i.e., |p, i) = |p, ef)) and satisfy the normalization condition 



(p / ,i|p,2) = (2 7 r)^(p-p / K 



'.i 



(3.15) 



In terms of these states the nonrelativistic expansion of the matrix element of the electro- 
magnetic current is (up to NNLO) 



\P',J\Jem\P,i) = e 



F c{q 2 )^ + ^2 F q(i 2 ) (w 



1 



n 



-<i 2 5ij 



2M d . 
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<p',j|J*Jp,i> = ^rr [W)fy(P + P? + F M (q 2 ) (#<* - tfq, 



2Mi 



a^ofo 2 ) ( W - 7737<l%j OP + P? 



(3.16) 



?o — | <^l| 2 , is the square of the four-momentum transfer, and n is 



where q = p' — p, q 2 

the number of space-time dimensions. (Note that we are using a different normalization of 

states to JL^,|(J.). The dimensionless form factors defined in eq. ( [3.16D are normalized such 

that 

F c (0) = l , 



2M, 



■F M (0) = Hm 



_F Q (0) = \x Q 



(3.17) 



where [Lm = 0.85741 ^— is the deuteron magnetic moment, and /ig = 0.2859 fm 2 is the 



(3.18) 



deuteron quadrupole moment. The charge radius of the deuteron J (r 2 .) is defined by 

Fc(q 2 ) = l + \(r 2 d )q 2 + ... . 



In effective range theory, the short-distance part of the deuteron wave function is only 
important for establishing the charge normalization condition, Fq{0) = 1. The prediction 
of effective range theory for the form factor Fc(q 2 ) follows from the Fourier transform of 
|^*>(r)| 2 , 



C 



(q 2 ) 



1 



4 7 



This yields a deuteron charge radius, 

„2\ER _ ' ' 



tan 



-1 



47 



ipd 



(3.19) 



(r 2 d Y 



S7 2 1 - 7p rf 
1 



*r 



2^2 



1 + ipd +7 p d + 



(3.20) 



which gives a numerical value of J (^) ER = 1.98 fm, very close to the quoted value of the 
matter radius of the deuteron [0,0 of r m = 1.967 ± 0.002 fm. Conventionally, the charge 
radius is obtained by combining the nucleon charge radius in quadrature with the matter 
radius, and agrees very well with the experimental value. 



In the EFT(/f"), the charge form factor Fc(q ) has an expansion in powers of Q, F ( 



c\ 



Fc (q 2 ) + Fq {q 2 ) + Fq (q 2 ) + ... where the superscript denotes the order of the contribution. 
The leading contribution to charge form factor F c {q 2 ) is calculated to be 



,(o) 



F&'tf 



4 7 , -1 
tan 



V=¥ 



47 



(3.21) 
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which reproduces the leading term in a 7p^ expansion of the form factor computed with 
effective range, eq. ( 3.19Q . At NLO, we find a contribution to the charge form-factor of 



*£V 



-Ipd 



47 



V=q- 



tan 



-1 



47 



(3.22) 



which reproduces the subleading term in the jpd expansion of the effective range result, 
eq. Q3.19Q . At NNLO, there are contributions from several types of graphs, including rela- 
tivistic corrections. We find 
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32 

47 
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: tan 
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47 



(3.23) 



where the first line of eq. ( |3.23| ) reproduces the subsubleading term in the effective range 
expansion, eq. ( p.!9| ), while the second term is the relativistic correction. The third term is 
the contribution from the isoscalar charge radius of the nucleon, J(r% ), which is measured 



to be J{r% t0 ) = 0.79 ± 0.01 fm. This measured value contains both relativistic corrections, 

0.18 fm (the Foldy term), and contributions from strong interactions. 



, r 7V,rel) 



AM% 



There is no reason to separate these two contributions, as (r% ) is the coefficient of the 
nucleon-photon charge radius operator in the single nucleon sector. 
The charge radius of the deuteron resulting from this form factor is 



(ri) 



2\EFT 



5/y* 



( r ^,o) + 377 l + 7Pd + 7 



2^2 



+ 



32M^ 



0.62 + 
4.30 fm 2 



2.33 



0.95 + 0.39 + 0.0014 



(3.24) 



Taking the square root of this value 
2.07 fm, which is within a few percent of the measured value of J (r^) = 



where the last term is the relativistic correction, 
gives y^)" 

2.1303 fm [49,51 -Q. The magnitude of the relativistic correction we have computed, 
( r drei) = +0.0014 fm , is the same as the contribution from the spin-orbit interaction (rela- 
tivistic effect) computed in f52f , but is of opposite sign. 

A comparison between eq. ( |3.24|) and eq. ( |3.20|) reveals that charge radius computed up 
to NNLO in EFT (ft) is the same as that computed with effective range theory up to very 
small relativistic corrections. However, at N 3 LO there is a contribution to the charge form 
factor and to the charge radius from a four-nucleon-one-photon operator of the form 

O ct = (N T P l N)\N T P i N)v 2 A . (3.25) 

Therefore at N 3 LO the prediction of EFT(yf) and effective range theory will differ. 
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C. Magnetic Form- Factor of the Deuteron 

The magnetic form-factor of the deuteron has been computed in the effective field theory 
with dynamical pions |I^| and is found to be dominated by the nucleon isoscalar magnetic 
moment. At NLO there are no mesonic corrections other than those that can be absorbed 
into the definition of the nucleon magnetic moment. However, there is a contribution from 
a four-nucleon-one-photon counterterm with an unknown coefficient. As the pions play 
no explicit role in the deuteron magnetic moment at NLO, one can immediately write an 
expression for the magnetic moment and form factor in EFT^) from the expressions given 

ra. 



m 



The Lagrange density describing the magnetic interactions of the nucleons is 



A 



B 



2M, 



N 



■rf (k + «iT 3 ) o ■ BiV 



(3.26) 



where k 



K, 



, + K n ) and K\ = h( K v ~ K n) are the isoscalar and isovector nucleon magnetic 



moments in nuclear magnetons, with 

Kr, = 2.79285 , 



-1.91304 



(3.27) 



The magnetic field is conventionally defined B = V x A. At NLO there are four-nucleon-one- 
photon operators that appear and can contribute to both the deuteron magnetic moment 
and the rate for np — > dj, 



C 



2,-B 



e % (N T P i N)\N T P 3 N)B i - e*L 2 ie ijk {N T P i N)\N T P j N)B k + h.c.] . (3.28 



At NLO the deuteron magnetic moment is 



fJ<M 



K 



p + K>n + L^ 



2Mjv7(/i - 7)' 



2M N \ ' ir 

Comparision with the measured value of hm gives 

%(m n ) = -0.149 fm 4 , 



(3.29) 



(3.30) 



at the renormalization scale \x = m*. The evolution of the ^(/-O operator as the renormal- 
ization scale is changed is determined by the RG equation 
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dfi 



0,-1 



The magnetic form factor has an expansion in powers of Q, -FAf(<? 2 ) 
f m(q 2 )+ •••• We find that 
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(3.31) 
Ffhf) + 

(3.32) 



and 
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FIG. 3. The form factor A{q 2 ) as a function of |q| = y 
the leading order prediction, the dotted curve corresponds to the next-to-leading order prediction, 
and the solid curve corresponds to the next-to-next-to-leading order prediction, in EFT(/f). 
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+ e*L 2 ^(/i- 7 ) 



(3.33) 



Despite the fact that ^L 2 is numerically small, it is important to realize that the effective 
range expansion for the deuteron magnetic moment is not formally valid beyond leading 



order (where jj,m 



[k p + K n )). At NLO, there is a counterterm that is allowed by the 



2Mjv *• p 
symmetries of the theory and it is expected to be of natural size, which in EFT(yf) is set by 

the pion mass. Given that there is a counterterm at NLO, we do not pursue this calculation 
to higher orders, even though it is straightforward. 

It is combinations of the electric, magnetic and quadrupole form factors that are mea- 
sured in elastic electron-deuteron scattering. The differential cross section for unpolarized 
elastic electron-deuteron scattering is given by 



da 
dfi 



da 
dfi 



Mott 



A(q 2 ) +B(q 2 ) tan 2 (? 



(3.34) 



where A and B are related to the form factors that appear in eq. (|3.16|) by |^9 

2 



A 



F 2 C + -^Fl + ^F 2 Q 



4 
—1 



B = ^(l + V )F 2 M , 



(3.35) 



with i] = — g 2 /4M|. In order to compare with data, we take our analytic results for the form 
factors and expand the expression eq. ( |3.35|) in powers of Q. At the order we are working, 
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FIG. 4. The form factor B(q 2 ) as a function o/|q| = \J —q 2 ■ The dashed curve corresponds to 
the leading order prediction, and the dotted curve corresponds to the next-to-leading order predic- 
tion, in EFT(/f). 



A is sensitive both the electric and magnetic form factors, while B depends only on the 
magnetic form factor. The predictions for A(q 2 ) and B(q 2 ) along with data are shown in 
fig. (H) and fig. (|4]), respectively. 



D. Electric Quadrupole Form-Factor of the Deuteron 

The quadrupole form factor is dominated by mixing between the S-wave and D-wave 
components of the deuteron due to the NN interaction. However, at subleading orders there 
are contributions from two-nucleon-one photon operators that do not contribute to nucleon- 
nucleon scattering (the same as for the deuteron magnetic moment). In the theory with 
pions, the leading contribution to the quadrupole form factor is from the exchange of a 



potential pion [Eq|, occurring at NLO in the power counting. The two-loop graphs give 
'(|k| 2 ) _ 3g 2 A M Nl fi . 1 



Ft ] > 
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dx 
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One finds the quadrupole moment at this order to be 






g 2 A M N (6-f 2 + 9 1 m n + 4ml] 
307r/ 2 (m T + 2 7 ) 3 



(3.36) 



(3.37) 



(3.38) 



Numerically, this is liq = li q = +0.40 fm 2 , about 30% larger than the measured value. 
This magnitude of deviation is expected based on the size of the expansion parameter. 



Subleading contributions, at NNLO in the theory with pions, have been computed in [54 



They depend upon the value of a coefficient that has not yet been determined from nucleon- 
nucleon scattering in the theory with pions, and so a numerical prediction at this order does 
not presently exist. Notice, that the expression in eq. (|3.38 ) is of order Q' 1 in the theory 
with dynamical pions, but is of order Q° in EFT(/r). 

In EFT(/f), the contribution from single and multi-pion exchange and from all other 
meson exchanges to the mixing between the S-wave and D-wave are captured in the ^Cq 
and Kj2 operators, up to NLO. However, at NLO there is also a contribution from a 
four-nucleon-one-photon operator, described by the Lagrange density 
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-e tl Q (N T P i N)^{N T P j N) (VV — V 2 <^') A 



(3.39) 



The operator connects two nucleons in initial and final S-wave states to a D-wave photon. 
As we have done with the other form factors, we expand the quadrupole form factor in 
, , t j. q r .... The leading order form factor is found to be 



powers of Q, F Q = F$ ] + F& + 
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At subleading order, we find a contribution to the quadrupole form factor of 
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(3.41) 



where we have written the form factor in terms of the three-momentum transfer and not the 
four-momentum transfer as we do not encounter relativistic effects at this order. 
At LO and NLO the deuteron quadrupole moment, //g, is found to be 
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Numerically, one finds that the quadrupole moment is given by the sum of fig = 0.273 fm 2 
and /i^ = (0.056 - 0.0165 t! Q ) fm 2 giving, at NLO, 



/i Q = (0.329 - - 0.0165 t Q ) fm 2 



(3.43) 



where ^Cq is measured in fm 5 . Setting ^Cq = 0, we find a quadrupole moment of 0.329 fm 2 , 
which is to be compared with the measured value of 0.286 fm 2 . In order to reproduce the 
measured value of the quadrupole moment ^Cq = +2.60 fm 5 , at jx = m^. 
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FIG. 5. The quadrupole form factor Fq(\]<.\ 2 ) as a function o/|k|. The dashed curve corresponds 
to the leading order prediction, and the dotted curve corresponds to the next-to-leading order pre- 
diction, in EFT(f). The dot-dashed curve corresponds to a calculation with the Bonn-B potential 
in the formulation of J5q/ . 



In contrast to the electric and magnetic form factors, the quadrupole form factor is 
found to deviate significantly from calculations in potential models (e.g. the model of |55| ) 
for |k|> m n , as shown in fig. (|J). However, within the region |k|< m^ where the EFT(f) 
is expected to perturbatively converge to the actual form factor we see deviations of less 
than ~ 10%, consistent with the size of the expansion parameter. It is clear that a NNLO 
calculation of this object is required in EFT(f). 



E. The Radiative Capture Process np — > dj 

The radiative capture process np — > d'y is a classic nuclear physics demonstration of 
the existence of meson exchange currents. In addition to the nucleon magnetic moment 
interactions which dominate this cross section, photons minimally coupled to pions have 
been estimated to contribute at the ~ 10% level in various schemes. The cross section for 
this process was computed in the EFT with KSW power counting in 



17 , where it was 



found that the pion exchange current contributions are ultra-violet divergent and require 
the presence of a counterterm at NLO (it would have been present at NLO even if the graphs 
were convergent). In EFT(/F), there are no dynamical mesons, and hence no meson exchange 
currents. For this process, and all the observables discussed in this paper, the effects of meson 
exchanges are reproduced order by order by local operators involving multiple nucleon fields. 
At NLO in EFT(/f) this corresponds to a single insertion of the *L\ operator, defined in 
eq. (p|). 

The amplitude for the radiative capture of extremely low momentum neutrons np — > d^ 
has contributions from both the ^So and 3 S\ NN channels. It can be written as 

iA(np -* drj) = e X N T r 2 a 2 [a- ■ k e(d)* ■ e( 7 )* - er ■ e( 7 )* k • e(d)*] N (3.44) 

+ ie Y e ijk e(dY* k j e( 7 ) fc * {N T T 2 r 3 a 2 N) , 
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FIG. 6. The Feynman diagrams giving the leading order contribution to np — > d"f in EFT(yf). 
The solid lines denote nucleons and the wavy lines denote photons. The light solid circles correspond 
to the nucleon magnetic moment coupling of the photon. The crossed circle represents an insertion 
of the deuteron interpolating field . 



where e = |e| is the magnitude of the electron charge, N is the doublet of nucleon spinors, 
5(7) is the polarization vector for the photon, e(d) is the polarization vector for the deuteron 
and k is the outgoing photon momentum. The term with coefficient X corresponds to 
capture from the 3 S\ channel while the term with coefficient Y corresponds to capture from 
the ^o channel. For convenience, we define dimensionless variables X and Y, by 



X 




Y 




(3.45) 



Both X and Y have the Q expansions, X = X^ + X^ + ..., and Y = Y<® +yW + ..., where 
a superscript denotes the order in the Q expansion. The capture cross section for very low 
momentum neutrons with speed |v| arising from eq. Q3.45 ) is 
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(3.46) 



where a is the fine-structure constant. 

At leading order in EFT^) the amplitudes receive contributions from the Feynman 
diagrams shown in fig. (^) and are 
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(3.47) 



where a^ So ' = —23.714 ± 0.013 fm, is the scattering length in the x 5o channel, and K\ is 
the isovector magnetic moment defined in eq. ( |3.26j ). At next-to-leading order, NLO, the 
contribution arising from the Feynman diagrams shown in fig. (0) and fig. (§) is found to be 

in 



F (1) = \^l Pd (l - 7a (1So) ) 
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FIG. 7. Graphs contributing to the amplitude for n + p —>■ d + 7 at subleading order due to 
insertions of the C2 operators. The solid lines denote nucleons and the wavy lines denote photons. 
The light solid circles correspond to the nucleon magnetic moment coupling of the photon. The 
solid square denotes a C<i operator. The crossed circle represents an insertion of the deuteron 
interpolating field . The last graph denotes the contribution from wavefunction renormalization. 
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(3.48) 



where rjj ° is the effective range in the ^So channel. We have not computed X^- 1 ' as it can 
only contribute at NNLO since X^ ' vanishes. The RG evolution of $L\ was discussed at 
length in [17], where it was made clear that its behavior is much different from *L,2, the 
counterterm for the deuteron magnetic moment. In the absence of pions we find 
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(3.49) 



in order that the cross section for NN — > NN'-f with the initial nucleons in the ^So channel 
and the final nucleons in the 3 S"i channel be independent of the renormalization scale at all 
energies. The analytic structure of the amplitude ensures that the capture cross section will 
be /^-independent, if NN — > NNj is \i- independent. 
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FIG. 8. Local counterterm contribution to the amplitude for n + p — > d + j at NLO. The solid 
lines denote nucleons and the wavy lines denote photons. The solid circle corresponds to an inser- 
tion of the "*L\ operator. The crossed circle represents an insertion of the deuteron interpolating 
field. 



The cross section for this process has been measured very precisely for an incident neutron 
speed of |v| = 2200 m/s to be a expt = 334.2 ±0.5 mb |6| 
at NLO, at this incident neutron speed, of 



In EFT (fi) we find a cross section 



a i = (287.1 + 6.51 *L X ) mb 



(3.50) 



where ^L\ is in units of fm and is renormalized /i = m^. Requiring ay to reproduce the 
measured cross section <j expt fixes ^L\ = 7.24 fm 4 . 

We see that even in the theory without dynamical pions, one is able to recover the cross 
section for radiative neutron capture at higher orders. It is clear that in this theory the four- 
nucleon-one-photon operators play a central role in reproducing the low energy observables. 



In the theory with pions, one can see by examining the contributing Feynman diagrams |]T7 



that in the limit that the momentum transferred to the photon is small the pion propagators 
can be replaced by l/m 2 , while keeping the derivative structure in the numerator. This 
contribution, as well as the contribution from all hadronic exchanges, is reproduced order 
by order in the momentum expansion by the contributions from local multi-nucleon-photon 
interactions. From the calculations in the theory with dynamical pions, the value of ^L\ is 
not saturated by pion exchange currents as these contributions are divergent, and require 

Therefore, estimates of ^L 1 based on meson exchanges 
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the presense of the L\ operator 
alone are model dependent. 

The effective range calculation of np — > d'-f was first performed by Bethe and Longmire 
32| and revisited by Noyes ||57|| . After correcting the typographical errors in the expression 



for a that appears in the Noyes article, the expressions in the two papers [32,57] are identical, 



a 



(ER) 



2na k\ 7 6 (a^) 2 qg> / 



v|M^(2 - 7 a( 35 i); 
which when expanded in powers of Q is 

a^ = ^[,l(l-W 1So) ) 2 

N 



+ 



7<2 



( 3 Si) 



7<2 



VSo) 



^rf 0) 



(3.51) 



+ \l { Pd ~ 4 So) ) «?(1 - 7 a (l5o) ) 2 + \l {pa + rf 0) ) «?(1 - ~fa^) + 



(3.52) 



From this expansion, one finds that 
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yER,(0) =Ki ^ _ 7a ( 1 5o)) xER,(0) = ^ (3>53) 

and 

F ER <« = i« l7 ( Pd - r£*>) (l - 7« (l5o) ) + ^17 (p, + rf 0) ) . (3.54) 

At LO in the EFT(^) expansion the cross sections agree, however, at NLO the expressions 
are very different. In addition to the counterterm that appears at this order in the EFT^), 
the contributions from the effective range parameters are found to disagree. In the EFT(/f) 
the local counterterm is renormalized by the short- distance behavior of graphs involving the 
C*2 operators and hence the effective range parameters in both channels. Given, this behavior 
it is no surprise that the effective range contributions differ between the two calculations. 
Given the explicit p-dependence in Y^ and the absence of such dependence in y ER '( 1 ) it 
is amusing to ask if there is a scale for which the expressions are identical, with ^L\ = 0. 
Indeed such a scale exits, 

7r ( l5 o) a ( 1 5„) _ pd 



which, by inserting the appropriate values, gives scale p ER ~ 144 MeV, coincidentally close 
to p = m n . 

IV. CONCLUSIONS 

An effective field theory without dynamical pions can describe all low-energy two-nucleon 
processes, and generally all multi-nucleon processes. We have shown in detail how EFT(/f) 
systematically reproduces nucleon-nucleon scattering and the interactions of the deuteron. 
Effective range theory is seen to emerge as the uncontrolled approximation to EFT(/f) where 
the multi-nucleon-external current local operators are neglected. For most observables, this 
is a small effect, while for the radiative capture process np — > d'-f this omission gives rise to 
deviations at the 10% level, which has been known for decades. Inclusion of the local four- 
nucleon-one- magnetic-photon interaction which enters at NLO in EFT(yf) introduces a free 
parameter ^Li, which can be fit to reproduce the measured cross section. As a consequence 
of this ultra-violet behavior, the NLO contribution to the cross section in EFT(f) with 
^L\ = is unrelated to that determined by effective range theory with the simple neglect 
of such contributions. While this does not satisfy our desire to make a parameter free 
prediction for this cross section beyond leading order, it does allow us to rigorously relate S- 
matrix elements for different processes based solely on the symmetries of QCD. The deuteron 
magnetic moment is another example of an observable where a counter term ^Li appears at 
NLO, but numerically this is found to be small. 

In working to higher orders in EFT(yf) one encounters relativistic effects, starting at 
NNLO. The contribution of relativistic effects on nucleon-nucleon scattering, the electric 
charge form factor of the deuteron and the electric polarizability of the deuteron were calcu- 
lated at leading order. Their impact upon the static properties of the deuteron are expected 
to be very small, of order ~ 7 2 /M^, and this was recovered. 
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We should take a moment to ask about the impact this work will have on our under- 
standing of potential model approaches to nucleon interactions. As we have emphasized, for 
most processes the effective range expansion provides a relatively good approximation. For 
some of the observables, such as the deuteron charge radius and the electric polarizability, 
this arises because counterterms only appear at very high orders in the momentum expan- 
sion, while for others, such as the np —>■ dj, and the deuteron quadrupole moment, there 
is a four-nucleon-one-photon counterterm that prevents effective range theory from doing 
better than ~ 10%. Such counterterms cannot be fixed by the nucleon-nucleon scattering 
amplitude no matter how precisely it is measured, but they can be determined from inelastic 
scattering processes. Thus, multi-nucleon processes involving external currents receive con- 
tributions from interactions beyond those related by gauge invariance to the nucleon-nucleon 
interaction. The form and magnitude of such interactions will, in general, depend explicitly 
upon the choice of the nucleon-nucleon interaction, as is made clear by the RG behavior of 
the ^Li counterterm. 
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